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Abstract 


This  repoit  presents  a  new  way  to  model  and  identify  single-link  flexible  arms  when  Qjulomb  fric¬ 
tion  is  present  in  the  joint  In  order  to  isolate  the  effects  of  this  nonlinearity,  the  arm  model  is  divided 
into  two  submodels:  motor  and  beam.  The  two  are  coupled  throu^  the  torque  at  the  base  of  the  beam. 

A  systematic  method  is  developed  to  obtain  the  dynamics  of  a  lumped-mass  flexible  beam,  and  some 
properties  are  obtained  from  this  model.  The  influence  of  the  payload  on  the  model  is  given  special 
attention.  This  way  of  modelling  is  later  extended  to  the  general  case  of  distributed-mass  flexible  arms. 

Identification  in  the  presence  of  Coulomb  friction  is  carried  out.  A  new  method  is  proposed  to  estimate 
this  friction  and  to  reconstrua  the  motor  frequency  response  from  experimental  results,  which  are  highly 
distorted  by  this  nonlinear  friction. 

Experimental  results  are  presented  and,  finally,  conclusions  are  drawn. 
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1.  Introduction 


1.1.  General  Introduction 

An  important  research  effort  has  been  canied  out  in  the  past  few  years  to  study  the  control  of  flex¬ 
ible  stnictufes  and,  in  particular,  flexible  aims.  Several  p^rs  have  appeared  on  this  topic  studying 
different  aspects:  [1,2]  are  examples  of  controlling  the  endpoint  position  using  state  space  techniques; 
[3-5]  used  different  schemes  of  adaptive  controllers  in  order  to  take  into  accoimt  changes  in  the  load. 
Classical  frequency  dmnain  techniques  to  control  a  flexible  ann  with  two  degrees  of  freedom  were  used  in 
[6].  However,  very  little  effort  has  been  devoted  to  the  control  of  flexible  anns  when  static  and  dynamic 
frictions  are  present  in  the  joints,  even  though  this  case  is  very  common  in  practice. 

Our  work  is  devoted  to  solving  or  at  least  minimizing  problems  due  to  nonlinear  friction,  which  may 
be  very  noticeable  in  lightweight  flexible  atmf .  or  in  flexible  arms  moving  at  low  speeds  and  accelerations. 
We  propose  a  completely  new  way  of  mode  Jing,  identifying  and  controlling  these  arms  that  produces 
simpler  controllers  than  the  other  methods.  Our  final  goal  is  to  build  and  control  a  3-joint  flexible  arm 
with  3  degrees  of  fineedom,  but  we  have  centered  our  work  of  the  first  year  on  the  single-link  case. 

The  results  of  our  research  during  this  first  year  are  presented  in  three  reports.  This  first  one  is  devoted 
to  modelling  and  identification.  The  second  me  proposes  a  general  control  scheme  and  compares  three 
specific  controllers  for  the  tip  position.  The  third  studies  a  control  scheme  to  deal  with  changing  payloads. 

IJL  Introduction  to  Part  I  - 

First,  we  consider  the  problem  of  modelling  a  special  class  of  single-link  flexible  arms:  lumped-mass 
flexible  arms.  They  consist  of  massless  flexible  structures  that  have  some  masses  concaitrated  at  certain 
points  of  the  beam  (see  Figure  1).  Only  translations  of  these  masses  produce  stresses  in  the  flexible 
structure;  their  rotations  do  not  generate  any  torque  in  the  beam.  So,  the  number  of  vibrational  modes  in 
the  structure  coincides  with  the  number  of  lumped  masses. 

Little  work  has  been  done  on  the  analysis  and  control  of  lumped-mass  flexible  stractures.  Book  [7] 
studied  the  case  of  two  rigid  masses  cotmected  by  a  chain  of  massless  beams  having  an  arbitrary  number 
of  rotation  joints.  Our  problem  differs  from  this  in  the  sense  that  our  structure  has  only  one  rotation  joint 
and  an  arbitrary  number  of  lumped  masses  along  the  structure. 

These  particular  straaures  are  studied  here  because: 

•  Their  dynamics  may  be  modelled  easily  compared  to  distributed-mass  flexible  arms. 

•  Interesting  properties  for  the  control  of  flexible  arms  are  deduced  from  their  dynamic  models. 

•  A  method  of  easily  controlling  these  arms  is  inferred  from  the  stracture  of  the  model. 

•  The  influence  of  changes  in  the  tip’s  mass  is  easily  characterized. 

•  Some  robots  and  robot  applications  can  be  reasonably  approximated  by  these  models. 

•  Given  a  distributed-mass  flexible  arm,  there  always  exists  a  mmcated  dynamic  model  which  is  of 
the  same  form  as  the  lumped-mass  flexible  atm  models,  and  which  reproduces  the  dynamics  of  the 
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measured  variables.  This  will  allow  us  to  generalize  the  above-mentioned  control  method  to  the 
case  of  distributed-mass  flexible  aims. 


The  influences  of  changes  in  the  carried  load  are  studied  in  this  report  Results  of  this  study  allow  us  to 
tune  the  controllers  to  the  particular  payload,  in  order  to  keep  the  response  of  the  system  to  changes  in 
the  position  reference  approximately  invariant 

Section  2  establishes  the  dynamic  model  of  these  beams,  some  properties  of  the  model  are  deduced 
in  Section  3  that  are  useflil  in  controlling  these  aims,  and  this  modelling  method  is  extended  in  Section 
4  to  the  case  of  distributed-mass  flexible  beams. 

A  new  method  of  identifying  the  dynamics  of  flexible  arms  in  the  presence  of  Coulomb  friction  (which 
is  a  strong  nonlinearity)  is  developed  in  Section  S.  It  also  provides  an  accurate  estimate  of  the  average 
value  of  the  Coulomb  friction  over  the  range  of  working  velocities.  This  method  allows  us  to  reconstruct 
the  motor  fiequency  response  from  experimental  results,  which  are  highly  distorted  by  this  nonlinear 
friction. 

In  Section  6  these  methods  are  a^iplied  to  a  class  of  single-link  li^tweight  flexible  arms  that  we  have 
designed  in  our  laboratory  for  experimental  purposes.  These  anus  are  very  lightweight;  they  can  perform 
quidt  movements,  but  the  fliction  torque,  which  is  comparable  to  the  coupling  torque  between  the  motor 
and  the  beam,  precludes  the  use  of  other  existing  identification  and  control  methods.  Experimental  results 
are  given  for  two  of  these  arms. 

Finally,  conclusions  are  presented  in  Section  7. 


2.  Lumped'Mass  Flexible  Arm  Modelling 

We  divide  the  model  of  our  flexible  atm  into  two  submodels:  the  first  one  describes  the  behavior  of 
the  motor;  the  second  one  describes  the  behavior  of  the  mechaiucal  structure  using  the  angle  of  the  motor 
as  its  inpuL  These  two  submodels  are  coupled  by  the  reaction  torque  of  the  beam  on  the  motor  (see  Figure 
2).  This  model  is  quite  different  flom  the  models  normally  used  in  the  control  of  flexible  arms,  which 
consider  the  applied  torque  as  the  input  to  the  beam  (Trackenbrodt  [8],  Low  [9]).  Our  model  has  some 
advantages  when  identifying  flexible  arms  with  friction  in  the  joints  [10]  (explained  in  Section  5),  and 
when  trying  to  compensate  for  friction  [11]  (explained  in  Part  11).  We  show  here  that  another  advantage 
of  our  model  is  that  it  allows  us  to  separate  the  dynamic-model  parameters  that  depend  on  the  geometry 
of  the  beam  from  the  parameters  that  depend  on  the  lumped  masses  of  the  beam.  Special  attention  is 
paid  to  this  issue  in  this  section  because  it  will  be  used  in  the  control  design. 


2.1.  Beam  modelling 

Qmsider  the  system  of  Figure  1.  It  represents  a  massless  flexible  beam  with  n  point  masses  distributed 
along  the  structure,  the  last  mass  being  located  at  the  tip  of  the  beam.  The  inertia  of  the  motor  is  included 
in  the  motor  submodel.  Let  m,-,  1  <  i  <  n  be  the  values  of  the  masses  and  li  the  distances  between 
consecutive  masses  t  -  1  and  i,  where  l\  is  the  distance  between  the  rotation  axis  of  the  motor  and  the 
first  mass.  Let  Li  be  the  distance  between  mass  rm  and  the  axis  of  the  motor.  We  assume  that  beam 
deflections  are  small  enough  so  that  the  distances  between  masses  rm  (measured  along  the  length  of  the 
beam)  are  equal  to  the  distances  between  the  masses’  projections  on  the  x-axis. 
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Wc  establish  two  coordinate  systems,  both  with  origins  at  the  motor  axis:  the  i  -  y  system  is  fixed  in 
space,  while  the  x—y  system  moves  with  the  motor  shaft  Thus,  the  y-cooidinate  of  a  point  represents  the 
deflection  of  the  point  on  the  beam  from  its  initially  straight  configuration,  corresponding  to  the  x-axis. 
We  denote  as  F(x)  and  T\x)  the  force  and  torque,  respectively,  at  this  poim:  i.e.  the  force  and  torque 
acting  on  the  beam  just  to  the  left  of  a  point  due  to  the  acticm  of  the  beam  just  to  the  right. 

An  external  force,  Fn,  and  torque,  7^,  applied  at  the  tip  of  the  beam,  are  also  considered  in  our  model. 
Fn  represents  the  component  of  the  resultant  applied  force  that  is  normal  simultaneously  to  the  beam  and 
to  the  joint  rotation  velocity  veaor.  Tn  is  the  componem  of  the  resultant  torque  normal  to  the  x  -  y  plane. 
These  two  quantitities  represent  the  interacdon  of  the  beam  with  the  enviromnent  and  may  be  produced 
in  several  ways,  a  few  of  which  are;  they  may  be  the  reaction  forces  when  following  the  contours  of 
a  surface,  or  the  reaction  of  the  next  joint  when  dealing  with  a  multilink  aim;  they  may  be  produced 
by  a  load;  they  can  even  represent  the  small  effects  of  friction  on  the  tip  (for  example,  in  flexible  aims 
mounted  on  an  air  table). 

Assuming  small  deflecdons,  a  massless  beam  can  be  described  by  the  stadc  deflecdon  relation: 


El- 


d^y 

db^ 


If  the  sdf&ess  £  •  /  is  constant  throughout  each  interval  of  the  beam,  the  deflecdon  in  the  interval 
[{'—  1,!]  is  given  by  a  third-order  polynomial: 


y-xix)  -  ai,o  +  K.M  •  (X  -  £,_t)  +  Ui,2  •  (X  -  Li_i)^  +  Ui,3  •  (X  -  (1) 

where  are  the  polynomial  coefficients,  different  in  each  interval,  and  Lq  =  0.  We  assume  in  the 
following  analysis  that  £  •  /  is  constant  throughout  the  beam. 


2J.J.  Geometric  equations 

Imposing  continuity  condidons  between  two  consecudve  intervals  up  to  the  second  derivadve  we  get 
3  •  (rt  —  1)  equadons: 


yXli)  =  y.+  l(0)  =>  Ui,o  +  «,-.l  •  /.•  +  «,-,2  •  +  Mi, 3  •  /i  =  Mi+1,0 


dyi*i 

dx 


(0)  =►  Mi.i  +  2  •  Mf,2  •  /i  +  3  •  Ui,3.^  =  Uf+1,1 


(2) 


§(/.)  =  ^(0)  «..2  +  3  •  Mi, 3  •  li  =  Mi., ,2 

for  the  joint  between  (/  -  1,/]  and  [i,/+  1],  1  <  i  <  n.  The  continuity  condidon  between  the  motor  and 
the  beam  gives 
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(3) 


yi(0)  =  ^(0)  =  0  “1.0  =  0,  ui,i  =  0, 

and  the  condition  of  the  applied  torque  at  the  tip  gives 

=  +  =  (4) 

The  set  of  geometrical  equations  is  completed  expressing  deflections  at  speciflc  points  i  in  terms  of 
these  polynomials: 


Uija  +  “.,1  •  /,  +  «i,2  ■  if +  Ui,3‘^i=  yiLi)-, 


I  <  i  <  n 


(5) 


2.1.2.  Dynamic  equations 

If  we  apply  Newton's  equations  to  the  n  masses  we  get: 

0<i<rt  (6) 

>«+l  .  ^ 

F(Li)  =  ~^mj .  ^^+Fn  0<i<n  (7) 

y=«+i 

where  ^Lj)  (=  y(Lj)  +  •  Lj)  is  the  arc  position  of  mass  j  with  respect  to  the  x-axis  {9m  is  the  angle  of 

the  motor).  Taking  into  account  that 


Ui,2  = 


T{Li-x) 

lE-r 


and  combining  (6)-(8),  we  get 


Ui,3  =  - 


6  El 


^  6  E  l 
dfi  mi 


■  (ui^  -  umj): 


where  is  defined  according  to  (8). 


1  <  r  <  n 


(8) 


(9) 


2.JJ.  Dynamic  model  of  the  beam 

Equations  (2)-(5)  allow  us  to  express  the  coefficients  Uij  as  linear  functions  of  y(L,)  and  r„.  m,.3  co¬ 
efficients  can  be  expressed  in  a  compact  form: 
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(  Ui2  '' 

U22 

( yiLi)  \ 
yiLi) 

II 

j 

,  (10) 

V  “-3  j 

y(U) 

\  T„ 

where  U  is  a  constant  matrix  that  depends  only  on  the  dimensions  of  the  beam  and  the  location  of  the 
masses. 

Denoting  as  9i  the  angle  between  the  x-axis  and  the  radial  line  from  the  origin  to  mass  i  (see  Figure 
1).  taking  into  account  that  Bi  =  and  substituting  (10)  into  (9),  we  get  n  linear  equations  of  the  form 


^  ^  +  +  \<i<n 


(11) 


(12) 


Equations  (11)  and  (12)  can  be  expressed  in  a  compact  form  as: 


A<-^=£/.t^e+e<?»]  +  7?r„+Q-F,, 


(13) 


where  M  =  diag(mi,m2,...  ,m„),0^  =  ,9n),A  is  an  n  x  n  constant  matrix,  and  B^V  and  Q 

are  constant  n  x  1  column  vectors.  In  particular,  ^  -  (0,0, ...,0,L~‘).  In  equation  (13),  A,B,V  and 
Q  depend  only  on  the  geometry  of  the  beam.  The  values  of  the  lumped  masses  only  influence  the  matrix 
M. 


Z2.  Motor  modelling 

The  dynamic  model  for  a  D.C.  motor  is  simple  and  can  be  written  as: 


K  i=J- 


S9„ 

df- 


+  V. 


dt 


+  Cj  +  CF, 


(14) 


where  K  is  the  electromechanical  constant  of  the  motor,  i  is  the  current,  J  is  the  polar  inertia  of  the  motor, 
V  is  the  dynamic  friction  coefficient,  C,  is  the  coupling  torque  between  motor  and  beam,  and  CF  is  the 
Coulomb  friction.  Taking  into  account  that  Ct  =  -7X0)  =  -E-l -I  -  u\2,  we  can  express  this  torque  as  a 
linear  function: 


Cl  =  'H  ■  Q  hn^\  •  9m  +  h«+2  • 


(15) 


where  X  =  (hi ,  /12, . . . ,  h„),  and  hi,  1  <  i  <  n  +  2,  are  parameters  that  do  not  depend  on  the  masses  along 
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the  beam. 


3.  Properties  of  the  Lumped-Mass  Model 

The  behavior  of  the  flexible  arm  as  described  by  equations  (13H1S)  is  represented  in  Figure  2.  We 
assume  that  F,  and  are  perturbations  to  the  system.  Then,  in  order  to  perform  the  analysis  and  design 
of  the  control  system,  we  make  =  7.  =  0.  Some  interesting  properties  of  the  lumped-mass  model  then 
follow: 

1.  From  (13),  using  Laplace  transforms  we  find  that  the  transfer  fimctions  9i(js)/9m(.s)  =  Gi(s)  for 
1  <  i  <  n  have  terms  only  of  the  form  <j<n>n  being  the  number  of  lumped  masses.  Then 
the  poles  and  zeros  always  come  from  factors  of  the  form  +  z),  z  being  a  complex  number  in 
general,  and  the  following  cases  are  possible: 

a)  z  €  92,z  >  C  two  conjugate  roots  on  the  imaginary  axis. 

b)  z  6  z  <  0  ^  two  roots  on  the  real  axis,  of  the  same  magnitude  but  opposite  sign. 

c)  z  €  C,  ^(z)  ^0  =>  sets  of  four  symmetrical  roots  with  respect  to  both  teal  and  imaginary  axes. 

The  mechanical  structure  is  marginally  stable  because  we  assume  there  is  no  flriction  along  the 
beam  (no  energy-dissipating  phenomena),  and  the  firictim  on  the  tip  is  treated  as  a  perturbation. 
So  the  poles  ate  always  of  case  a).  The  zeros,  however,  may  be  of  any  of  the  three  cases.  Case  b) 
appears  quite  often  and  produces  non-minimum  phase  systems.  Control  of  these  systems  presents 
s(xne  problems,  as  stated  in  [1]. 

2.  The  dilTercnce  between  the  orders  of  the  denominator  and  numerator  of  G,(s),  Vi,  is  at  least  two. 
This  is  a  consequence  of  flexible  arms.  If  the  numerator  and  denominator  were  of  the  same 
order,  it  would  mean  that  an  instantaneous  change  in  the  position  of  the  motor  would  produce  an 
instantaneous  change  in  the  position  of  the  ith  point  of  the  flexible  arm  (Initial  Value  Theorem  of 
Laplace  transforms;  see  Kuo  [12]).  This  is  not  possible,  however,  because  flexible  beams  need 
some  time  to  propagate  the  motion  along  the  stmcture. 

3.  Separation  Property,  It  was  mentioned  in  the  previous  section  that  the  influence  of  the  lumped 
masses  on  the  dynamics  of  the  atm  may  be  perfectly  separated  from  the  influence  of  the  geometrical 
dimensions  of  the  beam  (see  equations  (13)  and  (IS)).  This  happens  because  the  beam  is  modelled 
as  massless;  thus,  its  shape  is  given  by  static  deflection  equations  that  depend  only  on  the  position 
coordinates.  This  property  is  used  in  the  ctmtroller  design. 

4.  Zeros  Invariance  Property.  Zeros  of  the  Gn(s)  transfer  function  between  the  tip  position  and  the 
angle  of  the  motor  9„  remtun  constant  independent  of  the  payload  of  the  tip.  This  is  easily  proved 
taking  into  account  that: 

e  =  (M-^-E  l-A)-^  ■E  I  B  9„,  (16) 

and  that  the  inverse  of  a  square  non-singular  matrix  A  is  given  by  ([13]).  Then  the  nth 

row  of  (M  ■  ^  -  E  ■  I  ■  A)~^  is  independent  of  m„,  as  are  the  zeros  of  G„(s). 

5.  The  perturbations  F„  and  T„  affect  the  dynamic  behavior  of  all  the  masses  of  the  beam  structure 
linearly  (equation  (13)).  F„  does  not  appear  explicitly  in  (15);  this  is  because  the  force  applied  to 
the  tip  does  not  appear  in  the  geometric  equations. 


S 


4.  Extension  to  Distributed-Mass  Flexible  Arm  Modelling 


This  section  extends  the  proposed  modelling  method  to  the  case  of  flexible  aims  with  distributed  masses. 

The  following  lemma  is  proposed: 

Lemma:  Let  us  assume  that: 

1.  There  are  no  internal  energy-dissipating  phenomena  in  the  structure  of  a  distributed-mass  flexible 
arm.  Energy  dissipation  happens  only  in  the  joint  in  the  form  of  friction.  This  means  that  property 
1  of  Section  3  remains  true  here.  Property  2  is  also  true  because  the  motion  propagation  from  the 
base  to  the  tip  of  the  structure  is  not  instantaneous. 

2.  The  state  space  model  of  the  distributed-mass  flexible  arm  may  be  truncated,  keeping  in  the  model 
Just  the  first  few  representative  vibrational  modes,  and  neglecting  the  other  high-frequency  modes. 
This  truncation  is  always  performed  as  a  preliminary  step  in  the  design  of  controllers  for  flexible 
arms. 

3.  We  measure  the  position  at  n  points  of  the  beam.  We  choose  n  equal  to  the  number  of  vibrational 
modes  considered  in  the  truncated  model  of  the  beam,  N,,. 

Then  an  equivalent  lumped-mass  flexible  beam  model: 


^=A-9  +  B-0„  (17) 

C,  =  'rt-0  +  A^i  (18) 

can  be  d^ned  that  reproduces  the  dynamics  of  the  distributed-mass  flexible  beam  for  the  range  of 
frequencies  of  interest. 

This  model  represents  the  tnmcated  dynamics  of  the  distributed  mass  flexible  beam.  The  external 
perturbations  Tn  and  were  assumed  to  be  0  when  deriving  this  model.  Notice  that  now  A  and  B 
depend  on  the  distributed  mass  of  the  beam. 

The  following  procedure  is  proposed  to  obtain  this  model  from  experimental  measurements  on  a 
distributed-mass  system. 

Procedure. 

1.  Experimentally  obtain  the  frequency  response  between  the  angle  of  the  motor  and  the  angle  of  the 
tip. 

2.  Fit  a  transfer  function  ^,(,s)  that  relates  the  angle  of  the  tip  to  the  angle  of  the  motor  to  this 
frequency  data.  This  transfer  function  must  exhibit  properties  1  and  2  of  Section  3.  The  order  of 
the  denominator  divided  by  two  is  the  number  of  vibrational  modes  Nv  considered  in  the  model. 

3.  Experimentally  obtain  the  frequency  responses,  with  respect  to  the  angle  of  the  motor,  of  other 
(Nv  -  1)  ''ints  of  the  structure.  No  special  conditions  about  the  location  of  these  points  are 
required  a  priori. 
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4.  Identify  other  (N„  -  1)  transfer  functions  Gi(s)  from  this  data.  The  poles  are  common  to  the  whole 
system,  so  these  transfer  functions  will  have  the  same  denominator  as 

5.  The  poititm  of  the  model  described  by  equation  (17)  may  be  obtained  from  these  identified  transfer 

functions  (we  denote  &i  as  by  the  following  procedure: 

Let  &  €  be  the  column  vector  that  represents  the  measured  angles  0\,02,...  ,6Ny  at  the 

selected  points  of  the  structure: 


6i(s)  =  - ;  1  <  i  <  Nv  (19) 


where  5^  =  (  1 

^  s* 

... 

s  e 

Mi  =  (  Mi.O  Mi,l 

Mi, 2  • 

•  •  Mi/i*— 1  ^  > 

Mi  e 

S={Sq  6i  82  ...  8n^-x  l), 

Notice  that  Gi  is  of  this  form  because  of  properties  1  and  2  of  Section  3. 
Then: 


<f6i  (  0  (^i)  s 
dfi  ~  S-S 


•  e^(s). 


(20) 


If  itij  are  the  elements  of  ^4,  and  hi  the  elements  of  5  in  equation  (17),  then  substituting 
(19)-(20)  in  (17)  and  identifying  coefficients  in  the  numerator  we  get: 


The  coefficients  of  (17)  are  obtained  from: 


(21) 


^  XI 

(  H\  '' 

X2 

M2 

.4  = 

X3 

■ 

M3 

\XN.  / 

\  MW, 

where  rows  Xi  €  are  calculated  from: 

(  0  Hi)  -  ■6  =  (^xi  0  )  . 


(22) 


(23) 


6.  Expression  (18)  for  the  motor-beam  coupling  torque  Q  may  be  obtained  by  the  following  steps: 


(a)  Identification  of  the  linear  part  of  the  dynamics  of  the  motor,  transfer  function  Gm(s)  = 

(b)  Distortions  in  the  identification  because  of  the  Coulomb  friction  may  be  avoided  by  using  the 
procedure  described  in  the  next  section. 
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(c)  Denoting  Cds)  by  ^  •  ff„(s)  and  substituting  in  (14)  (after  making  CF  =  0)  we  get: 

JC  •  dc(s)  ■  dz\s)  =  (Js^  +  Vs)-dc(s)  +  flc(s). 

And  taking  into  account  that  de(s)  =  6 -S  (this  is  consequence  of  (15)).  we  detennine  J,  V  and 
ite(5).  K  is  obtained  from  catalogue. 

(d)  Equation  (18)  is  obtained  firom  polynomials  /le  and  de  by  following  an  algebraic  process  similar 
to  that  shown  in  step  5  of  the  procedure  for  obtaining  this  model  (page  10). 

The  only  constraint  to  consider  when  choosing  the  —  1  sensed  points  is  that  the  numerators  m  of  the 
identified  transfer  functions  at  the  ATy  points  should  be  linearly  independent  to  peimit  inversion  of  the 
matrix  that  tqrpears  in  (22).  This  means  that  very  close  points  carmot  be  used  for  sensing. 

A  consequence  of  this  section  is  that:  it  is  not  necessary  to  build  distributed-mass  flexible  arms  in 
order  to  test  control  schemes  for  flexible  arms.  Uonped-mass  flexible  arms  (which  may  be  easier  to  build, 
model  and  identify)  may  be  used  in  many  cases  as  prototypes  of  real  flexible  arms,  allowing  us  to  test 
control  laws  on  them  b^ore  implementation  in  real  arms. 


5.  Parameter  Identification  in  Cases  Involving  Coulomb  Friction  in  the 
Joint 

Several  methods  exist  that  identify  transfer  functions  of  dynamic  systems.  Because  of  the  characteristic 
of  flexible  arms  having  sharp  resonant  frequencies  in  the  normal  range  of  operation,  frequency  methods 
seem  to  be  better  suited  for  identifying  these  systems,,  and  have  been  used  by  many  researchers  [1,14]. 
However,  they  dealt  with  arms  without  friction  in  the  joints.  These  methods  can  be  easily  extended  to  the 
case  of  linear  dynamic  friction,  but  they  give  erroneous  results  when  the  Coulomb  friction  is  significant 
This  is  because  they  ate  based  on  linear  models,  but  the  Cbulomb  friction  is  nonlinear. 

In  this  section  we  develop  a  new  method  of  identifying  flexible  arms  that  extends  the  frequency 
identification  techniques  to  cases  involving  Coulomb  friction. 

The  dynamics  of  the  joint  arc  given  by  equations  (14)  and  (IS).  The  coupling  torque  C|(0  is  related  to 
the  angle  9m(t)  by  a  linear  differential  equation  (or  transfer  function)  obtained  by  as.<mming  Tn^Fn  =  0, 
using  Laplace  transforms,  and  substituting  Q  of  (13)  in  (IS).  The  term  CF(t)  is  related  to  the  velocity  of 
the  motor. 

Several  models  have  been  proposed  to  describe  the  friction  of  a  DC  motor,  taking  into  account  its 
linear  and  nonlinear  quantities  (a  recent  review  may  be  found  in  [IS]).  The  Coulomb  friction  is  the  most 
important  quantity  in  all  of  them.  This  is  especially  true  when  dealing  with  direct-drive  arms  because  the 
range  of  motor  speeds  is  relatively  low  and,  consequently,  the  dynamic  friction  is  also  low.  The  dynamic 
friction  may  be  modelled  as  a  first  approximation  by  a  linear  term  proportional  to  the  speed  of  the  motor, 
and  the  Coulomb  friction  CF  by  a  constant  whose  sign  changes  with  the  sign  of  the  velocity  of  the  motor. 
Thus  we  use  the  following  friction  model; 

Friction  torque  =  CF  +  V-  0„(t)  >  0 

Friction  torque  =  -CF  +  V-  0„(t)  <  0 

'Fhis  simple  model  is  a  good  description  of  the  friction  in  many  cases.  The  identification  method  that 
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we  propose  provides  an  averaged  value  of  all  the  parameters  (either  friction  or  linear  dynamics  parameters) 
of  the  aim  because  it  uses  the  spectral  characteristics  of  the  input  and  ou^iut  signals. 

The  proposed  identification  method  is  divided  into  three  stages  [10].  Hist,  a  hi^-gain  position  loop  is 
estabUsbed  around  the  motor  in  order  to  make  it  follow  the  reference  position  closely.  The  achievement 
of  this  allows  us  to  know  the  velocity  of  the  motor  and  hence  the  shape  of  the  temporal  evolution  of  the 
Qnilomb  fiicticn.  Then  get  the  magnitude  of  the  Coulomb  fiicdon  fiom  this  shape  and  from  the  spectral 
analysis  of  the  measured  signals.  Finally,  this  value  is  used  to  correct  the  experimoital  measurements  of 
the  fiequency  re^nse  of  the  motor,  which  are  distorted  by  this  nonlinear  Metion  term. 


5.1.  Control  oS  Motor  Position 

The  first  stage  of  our  proposed  identificatitMi  method  is  to  design  a  closed-loop  control  system  around 
the  position  of  the  motor  in  order  to  force  it  to  follow  a  specified  trajectory.  The  design  of  this  control 
system  is  straightforward,  and  a  normal  proportional-plus-derivative  (PD)  controller  may  be  used  for  this 
purpose.  If  we  apply  a  sinusoidal  reference  trajectory  to  the  closed-loop  system,  the  angular  velocity  of 
the  motor  will  be  r^roximately  sinusoidal.  The  motor  current  will  reflea  both  the  torque  needed  to  drive 
the  motor  and  beam  (i.e.  all  the  linear  terms)  and  the  torque  needed  to  overcome  the  Coulonb  friction, 
which  is  ideally  a  square  wave  in  phase  with  the  velocity  of  the  maor.  The  magnitude  of  the  Coulomb 
friction  remains  to  be  determined.  The  block  diagram  of  the  identification  setiq>  is  shown  in  Figure  3. 


SJZ.  Calculation  of  Magnitude  of  Coulomb  Friction 

Assume  that  the  closed-loop  position  control  system  for  the  motor  is  excited  with  a  sinusoidal  refer¬ 
ence  signal  of  frequency  uq.  If  the  system  were  linear,  we  would  find  that  the  current  and  the  angle  of  the 
motor  would  be  represented  by  sinusoidal  functiois  of  the  same  fiequency.  Because  of  the  nonlinearity 
present  in  the  system,  however,  the  current  will  not  be  sinusoidal.  Performing  a  spectral  analysis  of  the 
position  and  current  signals  ([16],  e.g.),  we  find: 

1.  The  motor  position  has  a  peak  at  fiequency  a;of  with  no  other  significant  peak. 

2.  The  motor  current  has  a  dominant  peak  at  frequency  u;o>  and  other  peaks  of  smaller  but  not  negli¬ 
gible  magnitude  at  frequencies  3  ■  u7o,S  •  ojq,!  •  b;o, . . .. 


An  explanation  of  the  second  phenomenon  concerning  current  may  be  given  in  terms  of  the  Fourier  series 
expansion  of  the  Coulomb  friction  torque,  represented  by  a  square  wave  of  frequency  a;o.  The  Fourier 
series  of  a  square  wave  can  be  written  as 


OO 

CF(f)  =  5~!a,  •  cos(i  ■  (4^  ■  t  +  if>), 

Csl 


(24) 


where 


A.= 


2CF 

i  •  TT 


•(l-(-l)‘), 


(25) 
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CF  is  the  magnitude  of  the  Coulomb  friction  torque,  and  is  the  i^ase  angle  of  each  hannonic.  which 
in  this  case  is  equal  to  0  Vi  (all  the  harmonics  arc  in  phase  with  the  fundamental  wave). 

Equation  (25)  indicates  that  Coulomb  friction  generates  odd  haimomcs  of  the  fundamental  frequency 
UQ.  Because  these  hannonics  appear  on  the  rig^t  hand-side  of  equation  (14)  and  the  angle  of  the  motor 
does  not  contain  these  hannonics,  the  current  i  which  appears  on  the  left  hand-side  of  equation  (14)  must 
contain  these  hannonics. 

Using  equation  (25).  we  can  now  determine  the  magnitude  of  CF,  based  on  the  study  of  the  third 
hannonic. 

If  the  system  is  excited  with  frequency  u;o.  then  in  equation  (14),  the  fundamental  component  of  the 
current  is  given  by 


ii(wo)  =  G-\ijQ) .  ejiuQ) 

iv  *  T 


(26) 


where  i'i(a;o)  is  the  polar  representatitm  of  the  first  harmonic  of  the  current  at  uq.  G(!jjq)  is  the  frequency 
characteristic  function  of  the  linear  pan  of  the  motor  model  at  frequency  (jq,  and  LB  is  the  phase  of  B. 
Notice  that  the  term  j  •  is  added  to  express  that  the  Coulomb  friction  leads  the  motor  position 

by  90  degrees:  i.e..  it  is  in  phase  with  the  velocity  of  the  motor. 

Also  with  the  system  driven  at  uq,  the  third  hannonic  of  the  current  is  given  by 

/3(3  •  wo)  =  G-‘(3  •  wo)  •  B^O  •  «D)  +y  •  4-^  •  (27) 

While  9m3(3  *  wo)  is  close  to  zero  for  arms  without  any  flexibility,  this  is  not  the  case  for  flexible 
arms.  When  3  ■  wq  is  close  to  one  of  the  natural  frequencies  of  the  beam,  the  position  control  system 
is  imable  to  completely  compensate  for  the  varying  beam  torque  due  to  the  oscillations  that  are  excited 
at  that  frequency.  The  control  system  does  eliminate  these  oscillations  sufficiently  for  us  to  consider 
that  the  motor  position  still  closely  follows  a  sinusoid  of  frequency  wq  and  thus  equations  (24)  and  (25) 
approximately  valid.  However,  the  small  high-frequmcy  ripple  that  Bm  now  presents  should  be  considered 
in  equation  (27)  because  the  factor  G~‘(3  -wq)  is  very  hi^  (GO' -w)  s  0  at  the  beam  resonant  frequencies). 

We  repeat  this  experiment  now  using  a  sinusoid  signal  of  frequency  3  •  wq  as  a  reference  for  the  motor 
position.  This  time  we  get  the  fundamental  component  of  the  current  as 

ii(3  •  wo)  =  G-‘(3  •  wo)  •  M3  •  wo)  +y  •  (28) 

iv  •  TT 

Equations  (27)  and  (28)  constitute  a  system  of  two  complex  equations  with  two  unknown  complex 
parameters,  G~*(3  -  wq)  and  CF.  In  theory,  the  value  CF  should  be  real,  but  because  of  small  errors  in  the 
measurements,  it  was  found  to  have  a  small  imaginary  component.  In  order  to  get  the  real  number  CF 
that  gives  the  best  approximation  for  equations  (27)  and  (28),  an  algorithm  that  minimizes  the  following 
cost  fimction  was  implemented: 
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Cost  =  U?(\  CF*ICF  1)  +  (ZC^)*, 


(29) 


where  CF*  is  defined  fiom  CF  as  follows: 

1.  For  a  given  value  of  the  Coulomb  fnction  CF,  obtain  G~*(3  •  mq)  from  (28). 

2.  For  the  above  G“*,  otxain  CF*  from  (27). 

In  other  words,  a  fimction  is  defined  between  CF  and  its  associated  cost  according  to  the  scheme: 


CF 


(28) 


(27) 


>  CF. 


(29) 


♦  CF 


♦  Cost  CF 


The  value  of  CF  that  minimizes  this  cost  fimction  was  obtained  using  a  simple  direct  search  method. 
Notice  that  (29)  represents  the  magnitude  of  the  ratio  ln(^),  which  is  a  complex  number.  This  cost 
fimction  was  chosen  to  give  equal  weight  to  errors  in  magnitude  and  phase. 

Different  values  of  CF  must  be  obtained  for  different  frequencies,  and  the  average  of  them  must  be 
taken  as  the  estimated  value  for  the  Coulomb  frictitm.  This  is  because  the  values  of  CF  may  differ 
considerably  depending  on  cjq  (we  had  variations  of  20%  around  the  mean  value  in  our  experiments). 


S3.  Correction  of  Frequency  Data 

Once  the  mean  value  of  the  Coulomb  friction  is  obtained,  the  frequency  response  of  the  linear  pan 
of  the  model  of  the  motor,  Giu),  can  be  expressed  (from  equation  (26))  as: 


G(u;)  = 


i((jj)  —  j  • 


(30) 


Equation  (30)  is  solved  using  the  experimental  data  obtained  previously  for  the  current  i(u)  and  the 
angle  of  the  motor  Oniu). 
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Therefore,  we  have  obtained  the  corrected  experimental  frequency  data  and  an  average  value  for  the 
Cbulomb  fiiction.  Now  a  transfer  iimction  may  be  fitted  to  this  data.  This  transfer  function  gives  the 
relationship  between  the  current  and  the  angle  of  the  motor. 

Also  obtained  from  the  experimental  was  the  transfer  function  between  the  angle  of  the  motor  and 
the  angle  of  the  tip.  This  way  of  representing  the  dynamics  of  the  system  is  quite  different  from  the  way 
used  in  the  other  approaches,  in  which  relations  between  current  and  angle  of  the  motor,  and  betweoi 
current  and  angle  of  the  tip  are  estaUished.  These  two  last  relatitms  are  influenced  by  the  Coulomb 
firiction,  and  the  measurements  taken  to  identify  them  are  therefore  distorted  by  this  nonlinear  quantity. 
Our  representation  has  the  advantage  that  identifying  the  transfer  function  between  the  angle  of  the  motor 
and  the  angle  of  the  tip  produces  a  very  clean  transfer  function  that  is  independent  of  motor  friction. 
Consequently,  nonlinear  terms  are  not  present,  allowing  the  use  of  conventional  flequency  identification 
methods  to  obtain  the  dynamics  of  the  beam. 


6.  Examples  and  Experimental  Results 


6.1.  Experimental  Setup 

The  meehnniral  system  consists  of  a  DC  motor,  a  slender  arm  attached  to  the  motor  hub,  and  a  mass  at 
the  end  of  the  atm  floating  on  an  air  table.  Hgure  4  shows  the  major  parts  of  the  system.  Optionally, 
masses  at  intermediate  locations  of  the  beam  may  be  added. 

Hgure  5  shows  the  arrangement  of  the  single-mass  beam  used  in  the  first  example.  In  this  case,  the 
arm  is  a  piece  of  music  wire  (7  inches  long  and  0.032  inches  in  diameter)  clamped  to  the  motor  hub. 
The  tip  mass  is  a  l/l6-inch  thick,  S  3l^inch  diameter  fiberglass  disk  attached  at  its  center  to  the  end  of 
the  beam  with  a  fieely  pivoting  pin  joint  The  disk  has  a  mass  of  54  grams  and  floats  on  the  horizontal 
air  table  with  minimal  fliction.  Because  the  mass  of  the  beam  is  small  compared  to  that  of  the  disk, 
and  because  the  pin  joint  prevents  generation  of  torque  at  the  end  of  the  beam,  this  mechanical  system 
behaves  practically  like  an  ideal,  undamped  mass-spring  system,  being  a  minimum  phase  system. 

In  the  second  example,  the  wire  is  replaced  by  a  longer  one,  and  two  masses  are  attached,  one  at  the 
middle,  and  the  other  at  the  end  of  the  wire.  The  secmid  system  has  the  characteristics  of  flexible  arms 
with  distributed  mass,  being  a  non-minimum  phase  system. 

An  Inland  direct-drive  motor  drives  the  atm.  The  amplifier  current  limit  is  set  at  4.12  amp,  which 
corresponds  to  a  9.0  Ib.inch  motor  torque.  The  Coulomb  friction  of  the  motor  is  about  .288  Ib.inch 
(corresponding  to  .132  amp)  and  has  a  significant  effect  on  the  control  when  the  torque  to  the  arm  is  low, 
as  with  our  very  slender  arms. 

Two  sensors  are  used  for  the  control  of  the  system.  A  11%-inch,  360~degree  potentiometer  measures 
the  angle  of  the  motor  shaft  A  Hamamatsu  tracking  camera  senses  the  x-y  position  of  an  irtfrared  LED 
mounted  on  the  tip  of  the  atm  in  the  case  of  the  single-mass  flexible  atm.  In  the  case  of  the  two-mass 
flexible  arm,  a  Sel^t  tracking  camera  is  used  that  simultaneously  senses  the  position  of  two  infrared 
LED’s  mounted  on  the  two  masses  of  the  arm. 

The  control  algorithm  was  implemented  on  an  Omnibyte  OB68K1A,  MC  68000-based  computer  with 
512  Kbytes  dynamic  RAM  and  a  10-MH  clock.  Because  the  control  computer  is  relatively  slow,  real-time 
computations  were  done  in  integer  or  short  integer  mode.  Analog  interfacing  was  provided  with  12-bit 
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A/D  and  D/A  boards.  The  sampling  period  used  to  control  both  arms  was  3  msec. 


6^  Single-Mass  Flexible  Arm 
62.J.  Modelling 

This  is  the  simplest  case.  Fbilowing  the  procedure  in  Section  2,  we  have: 

Geometric  equations: 

•  There  is  no  equation  (2). 

•  Equation  (3):  ui,o  =  ui,i  =  0. 

•  Equation  (4):  +  3  •  /i  •  =  Ttl- 

•  Equation  (5):  ui,o  +  h  ■  “i.i  +  /j  •  +  /f  •  =  )<Li) 

Dynamic  equations: 

•  Equation  (9):  /i  •  ^  ^  •  (m  j  +  jfc) 

Dynamic  model  of  the  beam: 

•  Equation  (10):  ai^  =  (  )• 

u 

•  Substituting  y(Li  )  =  l\-(6i-  in  the  above  equation,  and  combining  it  with  the  dynamic  equation, 
we  obtain  equation  (13): 

yfio  'X  1 

(31) 

Dynamic  model  of  the  motor 

•  We  use  (14)  as  it  is. 

•  Coupling  torque  (15): 

C,  =  -2EIuia  =  6E/-liUu-Ti  =  -^-^&i  +  ^-^-9„-t-\Ti  (32) 

<1  <1  2 

There  is  dynamic  friction  between  the  air  table  and  the  disk.  Therefore,  F\  =  -v-  u  being  the  friction 
coefficient.  Assuming  T\  =  0,  and  substituting  the  mechanical  parameters  £  =  30  •  10®  lblin.^,1  = 

=  0.121  lb.,  and  /i  =  7  in.  in  equations  (31)  and  (32),  we  get: 
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(33) 


gi(.y) _ 43J _ 

9„{s)  +  1.18  •  •  5  +  43.1 


C(  =  0.662  lbAn./rad.‘(0„  —  di) 


(34) 


622.  Identification 

Transfer  functions  of  the  motor  and  beam  were  obtained  using  the  method  described  in  Section  5.  Figures 
6  and  7  show  the  frequency  responses  of  motor  and  beam,  respectively.  The  transfer  function  fitted  to 
the  corrected  data  of  the  motor  was: 


dm(s)  394.94 .  (5^  +  0.06  -5  +  43 .75) 

Us)  ■  5  •  (5^  +  2.26  •  52  +  165.7  •  5  +  103.56) 

The  transfer  function  fitted  to  the  data  of  the  beam  was: 


52+0.06.5+43.75’ 

From  these  transfer  functions,  the  following  parameters  were  obtained: 
Parameters  of  equation  (14): 

/  =  0.005529  lb.in.se<P' 

V  =  0.01216  lb.in.frad.lsec. 

K  -  2.184  Ib.in.lamp.  (value  obtained  from  catalogue) 
CFIK  =  Q.\31  amp. 

Equation  (15)  gives: 


Ct  =  0.674  lb.in./rad.  ■  (0m  - 


Equation  (13)  gives: 


(35) 


(36) 


(37) 


0.12136  •  ^  =  -5.3095  •  0x  +  5.3095  •  0m  -  0.0072816  • 

dr-  dt 

where  the  last  term  on  the  right-hand  side  of  the  equation  represents  the  perturbation  because  of  the 
friction  between  the  disk  attached  to  the  arm  dp  and  the  air  table:  v  =  0.051  lb. sec.  j rad. 

Comparing  (33)-(34)  with  (36)-(37),  we  see  that  errors  between  theoretical  and  identified  parameters 
are  less  than  2%.  So,  the  mechanical  model  of  the  beam  may  be  accurately  obtained  from  the  theoretical 
parameters  obtained  by  using  the  equations  in  Section  2. 
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Parameteis  J  and  K  of  equation  (14)  can  be  obtained  from  a  catalogue,  but  V  and  CF  have  to  be 
experimentally  estimated.  Figure  8  shows  friction  characteristics  of  the  motor  obtained  by  using  the 
classical  method  of  plotting  the  motor  velocity  versus  the  current  needed  by  the  motor  to  nm  at  steady 
velocity.  Because  of  the  stiction  present  in  the  system,  the  values  of  the  currents  at  low  velocities  could 
not  be  obtained.  Using  a  linear  extrapolation  of  the  mqxrimental  data  given  in  Figure  8,  the  value  of  the 
current  requited  for  the  Giulomb  Motion  toique  was  found  to  be  0.161  amp  for  the  positive  velocity  and 
0.152  amp  for  the  negative  velocity.  These  values  are  differem  from  the  value  obtained  from  the  proposed 
spectral  method.  Some  motor  control  experiments,  which  will  be  described  in  Part  n,  showed  that  the 
value  estimated  from  our  method  was  more  accurate  than  these  values  obtained  Man  the  classical  method. 
This  is  because  Coulomb  and  dynamic  frictions  are  estimated  in  the  classical  method  by  extrapolating 
from  hi^-speed  measurements  that  are  outside  the  working  range  of  velocities  of  this  motor,  while  in 
our  method  these  values  are  estimated  from  low-velocity  measurements  inside  this  range. 


63.  Two-Mass  Flexible  Arm  Modelling 

We  develop  in  this  subsection  the  equations  of  the  beam  shown  in  Figure  9,  and  the  coupling  between 
motor  and  beam.  Identification  is  not  done  here  because:  1)  the  parameters  of  equation  (14),  including 
Coulomb  Metion,  remain  the  same  as  in  the  previous  example  because  here  we  use  the  same  motor, 
2)  in  the  previous  subsection  we  showed  that  differences  between  experimental  and  theoretical  beam 
trarLsfer  functions  are  small.  Thus,  the  beam  transfer  function  can  be  accurately  obtained  by  modelling 
the  mecharucal  system. 

Now  again  following  the  procedure  in  Section  2,  we  obtain: 

Geometric  equations: 

•  Equations  (2); 


Ml,0  +  Ul.l  •  A  +  Mu  •  +  Mu  •  /i  =  M2,0 


Mi,i  +  2  •  Uu  •  ^1  +  3  •  Uu  •  i\  =  M2,i 


Mu  +  3  •  «13  •  li  =  U22 

•  Equation  (3):  mi,o  =  mi,i  =  0. 

•  Equation  (4):  U22  +  3-  h-  M23  =  2^. 

•  Equations  (5); 

Mi,0  +  /i  ■  Mi,i  +  /?  •  Mu  +  fj  ■  Uu  =  y(Li) 
M2.0  +  h  ■  M2,i  +  4  •  M23  +  4  •  M23  =  y(L2) 
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Dynamic  equations; 


•  Equations  (9): 


Dynamic  model  of  the  beam: 


Equation  (10); 


,  Sh  e  .E  I  . 

h  •  -75-  - - (“u  -  “2,3) 

dr  mx 


„  ..  _6  E  I  ,  F2  , 

'2)  • -5^  - 


V  “23  J 


(  y(U) 
yiLz)  I, 
V  ^2 


(38) 


where 


U  = 


1 


3 


/,  ./2-(3-/i+4./2)  \  Ihfik 


•  Substituting  y(Li)  =  lx  •  (5i  -  9m),y{E2)  =  (/i  +  fe)  •  (^2  -  ^m)  in  the  above  equation,  and  combining 
it  with  the  dynamic  equations,  we  obtain  equation  (13): 


M 


E  I  [A  [  I +e-0J  +  ^  r2  +  Q  F2, 


(39) 


where 


/?-/2-(3-/t+4./2)  V  ^ 


-(/i  +  /2)V2 


2=1  ?  I, 

71^ 


/2-/2-(3  /i+4./2) 


-(6-/,+6-/2  +  2-^  +  2- jj)  (/i  +  /2)-(3  +  2-^) 
0/»-^3/2)/f  ,2/. 


and 


g-  6  [  /i  +  3./2  +  2.;J 

/?  /2  (3  /,+4./2)  V  -/?/(/! +/2) 
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Coupling  torque  (15); 


Cj  =  —2  •  E  ■  I  •  u\^  =  2-  E  ■  l\  •  ui,3  — 


2E-l-yi 


_  (9,  \  12-£./-(/i-h/2)  k 

\h  )  /i  •(3/1+4. /2)  "  (3/1+4-/2) 


(40) 


where 


n  = 


6EI-ili+lz) 

/2(3./,+4-/2) 


(  -(H-2.J;)  1 ), 


Substituting  the  mechanical  parameter  £  =  30  •  10^  lb/ in?',  I  =  in?,m\  =  mz  =  0.121 

lb.,  /i  =  /2  =  6  in.  in  equations  (39)  and  (40),  we  get: 


=  ( -»•' 
y  0.0 


0.12136  0 

0  0.12136 


-176.6032  110.377 

27.59425  -22.0754 


66.2262 

-5.51885 


-0.047619 
017857 


)■ 

)  ^  (  2  053  )  ^3  =  4. 106,  A4  =  -0. 1429. 


*  0.0833 

Assuming  Tz^Fz-  0,  we  obtain  from  (39)  the  following  transfer  functions: 


<?i(j)  _  545.7  •  (j2  +  106.10833) 

9„(s)  ^  i^  +  1637.1  -52  + 57903.3175 


his)  -45.465  •  -  1273.3) 

5«(j)  j^  + 1637.1 -5^  + 57903.3175 


(42) 


Theoretical  natural  frequencies  of  the  beam  were  obtained  from  the  poles  of  these  transfer  functions. 
They  are:  6.014  rad./ sec.  and  40.0116  rad./ sec..  These  frequencies  were  later  measured  experimentally, 
and  differences  between  them  and  these  theoredcai  values  were  smaller  than  5%.  Finally,  notice  that  the 
last  transfer  function  (equation  (42))  is  non-minimum  phase  exhibiting  a  positive  zero  at  the  theoretical 
position  of  35.683. 


7.  Summary  and  Conclusions 

This  report  has  discussed  the  modelling  and  identification  of  single-link  flexible  arms  with  lumped  masses. 
Most  flexible  arms  have  their  mass  distributed  along  all  the  structure,  but  there  are  some  applications 
in  which  a  model  using  lumped  masses  may  be  useful.  A  typical  example  is  the  case  of  a  lightweight 
flexible  arm  carrying  a  heavy  load.  Here,  only  one  vibrational  mode  is  normally  important,  allowing 
approximate  modelling  by  a  massless  beam  with  a  mass  attached  to  the  tip.  Another  reason  to  study 
flexible  arms  with  lumped  masses  is  that  some  properties  and  control  methods  for  these  systems  may  be 
easily  extended  to  the  distributed  mass  case. 


20 


A  new  way  of  modelling  these  systems  was  presented  in  Section  2.  In  order  to  deal  with  nonlinear 
Coulomb  friction,  two  sulxnodels  have  been  defined:  one  that  describes  the  behavior  of  the  motor  and 
includes  friction  nonlinearities,  and  another  linear  submodel  that  describes  the  mechanical  behavior  of  the 
beam.  Both  submodels  are  coupled  by  the  torque  at  the  base  of  the  beam.  In  order  to  make  the  model 
more  general,  perturbations  in  the  tip,  represented  by  a  torque  and  a  force,  have  been  considered.  These 
perturbations  allow  us  to  consider  the  effects  on  the  model  of  friction  in  the  tip,  changes  in  the  carried 
load,  or  reaction  forces  produced  by  the  next  joint  (in  the  case  of  a  multiple-link  flexible  arm). 

Some  properties  were  deduced  from  these  dynamical  models  in  Section  3.  The  most  interesting  one  is 
the  Separation  Property,  which  defines  the  influence  of  the  masses  on  the  beam.  It  permits  representing 
the  product  of  the  mass  i  and  its  angular  acceleration  by  a  linear  combination  of  the  deflections  of  the 
beam  at  the  points  where  there  are  masses.  This  linear  combination  depends  only  on  the  geometry  of  the 
beam. 

Modelling  was  generalized  in  Section  4  to  the  distributed-mass  flexible  arm  case.  Using  the  identifica¬ 
tion  procedure  described  in  Section  4,  the  dynamics  of  distributed-mass  flexible  arms  may  be  represented 
by  models  of  the  same  form  as  those  represented  by  equations  (13)-(1S).  Another  conclusion  drawn  from 
this  section  is  that  equivalent  lumped-mass  flexible  arms  may  be  used  as  prototypes  of  real  flexible  arms 
in  order  to  test  control  laws.  This  section  does  not  make  any  statement  about  the  influence  of  the  tip 
mass  on  the  dynamics  of  the  systertL  It  seems  that  distributed-mass  flexible  arms  may  be  considered  as 
the  limit  of  lumped-mass  flexible  arms  when  n  — »  oo.  But  arudytical  proof  of  this  has  not  been  provided 
yet  If  it  is  true,  then  properties  of  Section  3  (which  describe  the  influence  of  the  tip  mass  on  the  model) 
remain  valid  in  the  distributed-mass  case.  Further  work  is  needed  to  verify  this. 

Section  5  described  a  new  method  of  identifying  the  linear  part  of  the  model  of  a  motor  in  the  presence 
of  Coulomb  friction  (motor  submodel  of  a  flexible  atm,  described  by  equations  (14)-(15)).  In  this  method, 
the  Coulomb  friction  of  the  motor  is  obtained  from  frequency  data  instead  of  temporal  dau,  as  is  the 
case  in  other  existing  methods.  The  proposed  method  presents  two  advantages:  1)  Coulomb  friction  is 
estimated  for  the  normal  range  of  velocities  at  v^ch  the  arm  is  operating,  whereas  the  other  methods 
estimate  the  Coulomb  ffiction  by  extrapolating  from  high-speed  measurements;  2)  the  proposed  method 
gives  an  average  value  of  the  Coulomb  friction  over  a  range  of  velocities,  whereas  the  other  methods 
determine  the  friction  only  at  certain  speeds.  This  method  may  be  applied  to  any  DC  motor. 

Fmally,  some  modelling  and  experimental  identification  results  were  presented  in  Section  6.  Modelling 
has  been  applied  to  two  lumped-mass  flexible  arms  that  we  have  built  in  our  laboratory.  The  first  one  is 
a  minimum  phase  sin^e-mass  flexible  atm,  while  the  second  is  a  non-minimum  phase  two-mass  flexible 
atm.  The  identification  method  has  been  applied  to  the  single-mass  flexible  arm,  and  to  the  DC  motor, 
which  is  common  to  both  arms.  Experimental  results  agreed  with  the  theoretical  model. 
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Figure  3:  Identification  seti^. 


Figure  4:  Experimental  setup  of  the  single-mass  flexible  beam. 
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Figure  5:  System  nomenclature  for  the  single  -m&ss  flxin. 


Figure  6:  Frequency  characteristics  of  the  motor  submodel  (including  coupling  with  the  beam). 
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Figure  8:  Qiaiactcrization  of  friction  using  temporal  data. 
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Figure  9:  System  nomenclature  for  the  two>mass  aim. 


26 


References 


[1]  Cannon,  R.  R,  and  Schmitz,  E. 

Precise  Control  of  Flexible  Manipulators. 

Robotics  Research,  1985. 

[2]  Matsuno,  F..  Fiikushima,  S.  et  al. 

Feedback  Control  of  a  FlexiUe  Manipulator  with  a  Parallel  Drive  Mechanism. 

International  Journal  of  Robotics  Research.  Vol.  6,  No.  4,  Winter,  1987. 

[3]  Harahima,  F.  and  Ueshiba,  T. 

Adaptive  Control  of  Flexible  Arm  using  the  End-Point  Position  Sensing. 

Proceedings  Japan-Usa  Symposium  of  Flexible  Automation.  Osaka  (Japan),  July.  1986. 

[4]  Siciliano,  B.,  uan,  B.S.  and  Book,  WJ. 

Model  Reference  Adqitive  Control  of  a  One  Link  Flexible  Arm. 

25th  IEEE  Cortference  on  Decision  and  Control.  Athens  (Greece).  December,  1986. 

[5]  Rovner,  D.M.  and  Caimon,  R.H. 

Experiments  Towards  on-line  Identification  and  Control  of  a  Very  Flexible  One-Link  Manipulator. 
International  Journal  of  Robotics  Research.  Vol.  6,  No.  4,  Winter,  1987. 

[6]  Ower,  J.  C.  and  Van  de  Vegte,  J. 

Qassical  Control  Design  for  a  Flexible  Manipulator  Modeling  and  Control  System  Design. 

IEEE  Journal  of  Robotics  and  Automation.  Vol.  RA-3,  No.  5,  October,  1987. 

[71  Book,WJ. 

Analysis  of  Massless  Elastic  Chains  with  Servo  Controlled  Joints. 

ASME  Journal  of  Dynamic  Systems,  Measurements,  and  Control.  Vol.  101.  September,  1979. 

[8]  Trackenbrodt,  A. 

Dynamics  and  Control  Methods  for  Moving  Flexible  Structures  and  Their  Application  to  Industrial 
Robots. 

Proc.  of  the  5th  World  Congress  on  Theory  of  Machines  and  Mechanisms.  Montreal  (Canada),  1979. 

[9]  Low,  K.  R 

A  Systematic  Formulation  of  Dynamic  Equations  for  Robot  Manipulators  with  Elastic  Links. 
Journal  of  Robotic  Systems.  Vol.  4(3),  1987. 

[10]  Feliu,  V.,  Rattan,  K.S.,  and  Brown,  RB. 

Model  Identification  of  a  Single-Link  Flexible  Manipulator  in  the  Presence  of  Friction. 

I9th  ISA  Annual  Modelling  and  Simulation  Conference.  Pittsburg  (USA),  May,  1988. 

[11]  Rattan,  K.S.,  Feliu,  V.,  and  Brown,  RB. 

A  Robust  Control  Scheme  for  Flexible  Anns  with  Friction  in  the  Joints. 

1988  NASA  -  Air  Force  Workshop  on  Space,  Operation,  Automation  and  Robotics.  Dayton  (USA), 
July.  1988. 


27 


[12]  Kuo.  B.C 

Automatic  Control  Systems.  Pientice-Hall,  1982. 

[13]  Lancaster,  P.,  and  Tismenetsky.  M. 

The  Theory  of  Matrices.  Second  Edition  with  Applications.  Academic  Press,  1985. 

[14]  Hastings,  G.G.,  and  Book,  W.J. 

Vmificatioo  of  a  Linear  Dynamic  Model  for  Flexible  Rdx)tic  Manipulators. 

Proceedings,  1986  International  Conference  on  Robotics  and  Automation.  San  Francisco  (USA), 
April.  1986. 

[15]  Canudas,  C,  Astrom  KJ.,  and  Braun  K. 

Ada^nive  Friction  Compensation  in  DC-Motor  Drives. 

IEEE  Journal  of  Robotics  and  Automation.  Vol.  RA-3,  n-6,  December,  1987. 

[16]  Oppenheim,  A.  V..  and  Schafer.  R.  W. 

Digital  Signal  Processing.  Prentice-Hall,  1975. 


28 


